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ABSTRACT 

A study of the distribution of OH gas in the central region of the nearby active 
starburst galaxy M82 has confirmed two previously known bright masers and revealed 
several new main line masers. Three of these are seen only at 1665 MHz, one is detected 
only at 1667 MHz, while the rest are detected in both lines. Observations covering both 
the 1665 and 1667 MHz lines, conducted with both the Very Large Array (VLA) and 
the Multi-Element Radio Linked Interferometer Network (MERLIN), have been used 
to accurately measure the positions and velocities of these features. This has allowed 
a comparison with catalogued continuum features in the starburst such as Hn regions 
and supernova remnants, as well as known water and satellite line OH masers. Most of 
the main line masers appear to be associated with known Hn regions although the two 
detected only at 1665 MHz are seen along the same line of sight as known supernova 
remnants. 
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1 INTRODUCTION 

Astronomical masers are useful probes of the interstellar 
medium in star forming regions, both in the Milky Way 
and external galaxies. Masers usually occur under specific 
conditions found in star forming regions so the presence of a 
maser immediately provides information on the state of the 
ISM and acts as an indicator for ongoing star formation. 

Studies of star forming regions within the Milky Way 
have the advantage that the masers can be studied at high 
linear resolution. In our own Galaxy, OH masers have typ- 
ical luminosities of 2 x 10 21 Watts and sizes of 2 to 20 AU 
ijMinier et al.|[2002l ). At the other end of the scale are the 
so-called megamaser sources in active galaxies, sources a fac- 
tor of 10 7 brighter than ordinary Galactic masers, such as 
that discovered in Arp220 (IC4553; lB~aan et al.lll982l ). In be- 
tween these extremes of the pheno menon are masers such a s 
those detected previously in M82 |Weliachew et all 1 19841) . 
These have observed intensities of ~10 25 Watts and are 
some t imes referred to as "kil omasers" (e.g. lHagiwara et al.l 
l200ll ; iHenkel fc Wilson] 1 1990t ). Studies of kilomaser sources 
in individual galaxies, such as M82, provides a sample of 
masers which are all at an approximately equal distance and 
hence can be observed with the same linear resolution. 

This paper describes an investigation into the main 
line OH maser population in the nearby starburst galaxy 
M82 using data from both the Very Large Array (VLA) 
and Multi-Element Radio Linked Interferometer Network 
(MERLIN). The maser population is investigated and their 
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Table 1. Summary of the observations of OH masers in M82. 



positions and velocity distribution are compared to other 
known objects and the distribution of other gases. 



2 MASERS IN THE M82 STARBURST 

M82 is one of the closest, and therefore best studied, star- 
burst galaxies and radio observations of this galaxy are nu- 
merous. Extensive studies of the distribut ion and dynamics 
of ne utral hydrogen have been made (e.g. IWills et al.ll2000l . 
120021 ) and these have been c ompared with the molecular gas 
as traced by CO emission (|Shen fc Lolll995l ). The angular 
resolution of the CO observations is, however, comparatively 
low. Observations of the transitions of the OH molecule can 
be made with the same instruments at similar resolution to 
the previous Hi studies and so provide a better comparison 
of atomic and molecular gas. 

An OH maser was first detected in M82 us ing the Ef- 
felsbe rg telescope by Nguvcn- Q-Rieu et al.l (|l976l ). Using the 
VLA, IWcliachc w et al. I (|1984D detected two main line maser 
regions in the galaxy and noted that, although both spots 



2 Argo et al 



are an order of magnitude brighter than the most luminous 
Galactic maser, the emission from those in M82 may be com- 
ing from a number of smaller masers within each spot. The 
resolution of these observations however was not sufficient to 
detect any substructure or extension in either masing region. 

Subsequent observations have resulted in further de- 
tections of masers at other frequencies. A total of six OH 
features wer e detected using the 1612 and 1720 MHz satel- 
lite lines bv ISeaauist et al.l 1 1997h . Four features were seen 
in emission at 1720 MHz, two were seen in absorption at 
the same frequency, while only one feature was detected 
in emission at 1612 MHz at the same position as one of 
the 1720 MHz absorption featur es. H2O masers have also 
been detected i n M82 at 22 GHz (|Baudrv fc Brouillet|[l996l: 
lHariwaral[2005l ). 

Recently MERLIN and the VLA have been used to 
probe the OH absorption across the central starburst region 
in M82. As well as deep absorption features, eleven main 
line OH masers have been detected to a limit of 3cr, nine of 
which are new detections. 

Unless otherwise stated, all positions are given in J2000 
coordinates in the format bKbb corresponding to 

09 h 55 m aa!aaa and +69 °40'bb'.'bb respectively . At the dis- 
tance of M82 (3.2 Mpc. TFreedman et al.lll994r ), one arcsec- 
ond corresponds to a linear distance of 15.5parsecs. 
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53.11+48.0 


53! 094 


47'.' 88 


5.57 


<2.6 


<0.47 


50.97+45.3 


50.924 


45.37 


8.35 


7.83 


0.94 


50.35+44.1 


50.350 


44.13 


7.01 


29.9 


4.27 


50.02+45.8 


50.017 


45.83 


4.30 


5.09 


1.18 


49.71+44.4 


49.689 


44.17 


<2.6 


5.29 


>2.03 


48.45+42.1 


48.360 


41.93 


4.15 


6.94 


1.67 



Table 2. Definite maser detections made with MERLIN in 1995. 
The 3<t noise per channel was 2.6mjy/beam. The IDs are con- 
structed from the J2000 positions measured in the VLA 2002 
data-set. The flux densities are given in mJy and the listed ratios 
are for Si667:Si665- 
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RA 


Dec 


Sl665 
(mJy) 


Sl667 
(mJy) 


Ratio 


50.97+45.3 
50.35+44.1 


50?924 
50! 384 


45"37 
44'.' 24 


<18 
<18 


23.0 
57.3 


<1.35 
<3.36 



Table 3. Definite maser detections made with MERLIN in 1997. 
The flux densities are given in mJy and the ratio is Si667:Si665- 
The 3<7 noise was 17mjy/beam in the 1665 MHz observation and 
18mjy/beam in the 1667 MHz observation. 



3 OBSERVATIONS 

Data from the MERLIN archives, as well as more recent 
VLA observations, have been used to investigate the maser 
population in M82. A summary of the observations discussed 
here is given in Table fT] and further details are given in the 
following Section. 



in column one, is constructed from the VLA 2002 position 
where known. The exception to this is 50.02+45.8 whose ID 
is determined from the MERLIN positions listed in Table 
[5] as it was not detected in the 2002 data. Note that the 
spectra are shown smoothed with a Gaussian of width two 
channels, while the peak flux densities listed in Table [2] are 
measured from the unsmoothed data. 



3.1 MERLIN archive 

3.1.1 1995 observations 

MERLIN observations of OH in M82 were carried out on 
26th November 1995 using a total bandwidth of 7.75 MHz 
over 64 channels giving a velocity resolution of ~22kms _1 at 
a wavelength of 18 cm. Eight of the MERLIN telescopes were 
included in the array for these observations: Lovell, Mk 2, 
Tabley, Darnhall, Knockin, Defford, Cambridge and Wardle 
(the Wardle antenna has since been decommissioned). This 
provided a maximum baseline of 217 km and hence a reso- 
lution of (f.'2 at a wavelength of 18 cm. Gaussian fitting of 
components using the AIPS task jmfit resulted in a position 
uncertainty, relative to the phase calibrator, of ±0'.'02. 

This was a preliminary observation intended to detect 
the two known masers, look for additional sources of OH 
emission, and also to probe absorption by OH in the cen- 
tral starburst in M82 to provide a comparison with that 
detected in CO emission (for example). Due to the u-v cov- 
erage of MERLIN, however, extended absorption features 
are resolved out and very little absorption was actually de- 
tected in these observations. Several emission features were 
detected, however. Figure [1] shows spectra of the masers de- 
tected in this observation. Table [2] lists the masers found in 
1995 MERLIN data together with their positions and flux 
densities in this data-set. The ID label for each maser, given 



3.1.2 1997 observations 

These observations were carried out in May 1997 by switch- 
ing between 1612, 1665, 1667, and 1720 MHz using seven 
telescopes of the MERLIN array: Defford, Cambridge, 
Knockin, Darnhall, Mk 2, Lovell and Tabley. The correla- 
tor was configured to record 256 frequency channels 7.7 kHz 
in width, giving a velocity resolution of 1.4 km s~ . Only the 
1665 and 1667 MHz observations are discussed here. Figure 
[5] shows spectra of the two masers detected at 1667 MHz, 
illustrating the narrower frequency channel widths in this 
observation. Gaussian fitting components in these date re- 
sulted in aposition uncertainty of ±0'.'03. 

As these observations were carried out using narrower 
frequency channels than the 1995 observations there is a 
reduction in sensitivity. The 3er noise level per channel in 
these observations was 16.9mJy/beam at 1665 MHz and 
18.2 mjy/beam at 1667 MHz, the difference arising from the 
longer observation at 1665 MHz. 

As stated previously, the u-v coverage of MERLIN is not 
ideal for detecting extended absorption and most of it was 
resolved out in these observations. For detecting compact 
objects such as masers, however, MERLIN's resolution is an 
advantage as the observations are not greatly affected by 
deep absorption features. 

The detections made in the 1667 MHz data-set are listed 
in Table [3] No detections were made in the 1665 MHz data- 
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Figure 1. Spectra of the maser features detected in the 1995 observations with MERLIN. All are detected at >3<r in at least one line. 
The labels in the top right of each panel correspond to the IDs given in TablefS] The 3<r noise level in each channel is 2.59 mjy/beam and 
the spectra shown here have been smoothed in frequency using a Gaussian of width two channels. The ai-axis velocity scale is calculated 
relative to the rest frequency of the 1665 MHz line. 



set. In most cases in the other data-sets examined, the 
masers are significantly weaker in this line compared to the 
1667 MHz line. 



3.2 VLA 2002 

Observations were carried out in April and May 2002 us- 
ing the VLA in A configuration. At 1.6 GHz this gives an 



angular resolution of ~ l'/4. The data were taken in two 
observing runs, one of which was carried out in April, the 
other in May 2002. The data were recorded in 64 channels 
over a bandwidth of 6.25 MHz such that both main line 
maser frequencies were within the observing band over the 
extent of the velocity range within the disk of the galaxy. 
This setup resulted in a velocity resolution of 18kms _1 . 
3C286 (1331+305) was used for flux calibration, 0319+415 
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Figure 2. Spectra of the two maser features detected at 1667 MHz in the 1997 observations with MERLIN. The 3cr noise level in each 
channel was 18.2 mjy/beam. The data are shown smoothed in frequency using a Gaussian of width two channels. The x-axis velocity 
scale is calculated relative to the rest frequency of the 1667 MHz line. 



was used for bandpass calibration, and the phase calibra- 
tor was 0954+745. Calibration of the data was performed 
using standard techniques for the VLA. In order to isolate 
the line component, the radio continuum emission was sub- 
tracted using a set of line free channels near the start and 
end of the cube where the channel frequency fell outside that 
corresponding to the velocity range for OH in M82. 

Some of the maser features are bright enough to cause 
ringing in the 2002 VLA data. This effect can be largely 
removed by smoothing the data, although this results in 
adjacent channels no longer being independent. This was 
carried out using the AIPS task xsmth to apply a Gaussian 
function with a width of two channels. This also has the 
effect of reducing the measured peak flux since the emission 
is effectively smeared over several channels, hence the peak 
fluxes listed in Table U were measured in the unsmoothed 
data. 

The OH absorption across the central starburst in M82 
is well displayed in these observations, but the spectral res- 
olution of the data-set is not ideal for maser studies. De- 
spite this, several maser features were detected in emission 
in these data. Table [3] lists the positions of the maser fea- 
tures in this data-set which are >3cr, while Figure [3] shows 
a plot of the maser positions on top of a contour plot of the 
OH absorption across the M82 starburst. Gaussian fitting 
resulted in a position uncertainty of ±0'.'2. Figure 3] shows 
spectra of each maser feature listed in Table [4] 

The criteria used for classifying a signal as a maser was a 
detection of at least 3<r in two or more consecutive channels. 
This excludes some features which only appear in a single 
channel. Due to the wide channels of this observation it is 
possible that some masers are missed altogether. This could 
be because the emission is narrow but happens to fall at 
a channel edge, becoming spread over a couple of channels 
and hence too close to the noise level to be a firm detection. 

Another major factor affecting the detection of a maser 
signal in these data is the depth of the absorption. As can 
be seen in Figure [4] some of the masers are coincident with 
deep absorption features and are only detected because they 
are very bright. Weaker masers whose brightness is less than 
the magnitude of the depth of the absorption at the same 
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Figure 3. The location of the masers listed in Table [4] superimposed on a contour plot of the OH absorption in the 2002 VLA data cube 
which has been continuum subtracted and summed in frequency using the AIPS task SQASH. Contours are from -10 to +10 (in steps of 
1) X 0.2 mjy beam -1 . Solid contours are positive, dashed contours are negative. The diameter of each symbol is that of the restoring 
beam, l'/3. 



ID 


R.A. 


Dec. 


Sl665 


Sl667 


Ratio (1667:1665) 




(J2000) 


(J2000) 


(mjy) 


(mjy) 




53.62+50.1 


53?647 


50" 09 


2.44 mjy 


2.49 mjy 


1.02 


53.11+48.0 


53.105 


47.97 


1.86 


<0.77 


<0.41 


52.71+45.8 


52.707 


45.75 


5.02 


<0.77 


<0.15 


51.87+48.3 


51.868 


48.31 


1.71 


<0.77 


<0.45 


51.23+44.5 


51.230 


44.54 


4.85 


5.93 


1.22 


50.97+45.3 


50.974 


45.27 


13.0 


18.5 


1.42 


50.51+45.6 


50.513 


45.57 


1.85 


1.28 


0.692 


50.36+44.2 


50.360 


44.18 


10.7 


50.5 


4.72 


49.71+44.4 


49.707 


44.37 


<0.77 


2.34 


>3.04 


48.45+42.1 


48.445 


42.13 


2.64 


5.97 


2.26 



Table 4. Definite maser detections identified from the VLA 2002 data. Limits on the measured flux densities are given as 3cr. The maser 
IDs here are constructed from the positions measured in this data-set and are used throughout this paper. 



6 Argo et al 

position are therefore difficult to detect in these observa- 
tions. 



53.62+50.1 is coincident with an Hn region (44.93+63.9) 
and was detected in both main lines in all but the 1997 
data-set. 

53.11+48.0 is coincident with a supernova remnant 
(44.43+61.9) and, while it was seen at 1665 MHz in all but 
the 1997 data-set, it was not detected at 1667 MHz in any 
observation. 

52.71+45.8 is coincident with both a supernova remnant 
(44.01+59.6) and a satellite line maser (feature 4 from 
ISeaauist etafll 19971 ) . Like 53.11+48.0, this maser is not de- 
tected at 1667 MHz in any observation. It is detected at 
1665 MHz in the 2002 data-set. 

51.87+48.3 is apparently not coincident with a known con- 
tinuum feature or either of the other types of maser shown 
in Figure [7] This source was only detected in the 2002 data 
and was one of the weaker detections. 

51.23+44.5 is located within an arcsecond of an Hn region 
(42.48+ 58.4), and close to two wat er maser features (SI and 
P3 from iBaudrv fc BrouilletJ Il996l ) and two further Hn re- 
gions. Despite not being detected in either the 1995 or 1997 
observations, it is seen in both lines in 2002. 
50.97+45.3 appears coincident with both an Hn re- 
gion (42.21+59.2) and a water maser feature (S2 from 
IBaudrv fc Brouilletl Il996h . In the 1997 and 2002 observa- 
tions this source is brightest at 1667 MHz, however in 1995 
it was brighter at 1665 MHz. 

50.51+45.6 is not clearly associated with another feature 
in M82 and is brighter at 1665 MHz than 1667 MHz. 
50.36+44.2 is closely associated with an Hn region 
(41.6 4+57.9) and within a beam width of a water maser (S3 
from IBaudrv fc Brouilletl Il996l ). This source is the bright- 
est maser feature in the M82 starburst and is significantly 
brighter (by a factor of more than four) at 1667 MHz than 
it is at 1665 MHz, the most extreme line ratio out of all the 
detections. 

50.02+45.8 is located between a supernova remnant 

(41.29+59.7) and a possible water maser (P2 from 

IBaudrv fc Brouilletl [1999 ) and is about a beam width from 
each. It is only detected in the 1995 data-set where it is 
brightest at 1667 MHz. 

49.71+44.4 i s coincident with a possib le water maser fea- 
ture (PI from Ba udrv fc Brouilletl Il996l ) and also within a 
beam width of two Hn regions (40.95+58.8 and 40.96+57.9) . 
It is brighter at 1667 MHz in every data-set in which it is 
detected, not appearing above the noise in any data-set at 
1665 MHz. This maser appears to consist of two closely- 
spaced velocity components. 

48.45+42.1 is also within a beam width of an Hn re- 
gion (39.68+55.6), close t o a water maser feature (S4 from 
IBaudrv fc Brouilletl Il996l ) and is consistently brighter at 
1667 MHz. 



4 NOTES ON INDIVIDUAL SOURCES 

A brief description is given here for each source. A more 
detailed discussion of the properties of the maser popula- 
tion is given in the following section. Note that, while the 
masers are identified here according to their J2000 positions, 
the continuum feature s to which they are compared (from 
iMcDonald et al.ll2002j)_ are identified using the B1 950 coor- 
dinate convention of Kronberg fc Wilkinson! (1 19751 ). 



5 DISCUSSION 
5.1 Line ratios 

For an optically thin cloud in local thermal equilibrium, the 
ratio of emission for the ground states (1667:1665:1612:1720) 
is 9:5:1:1, so Si665>Si667 is impossible without excitation 
such as that provided by a maser pumping mechanism. Ta- 
ble [5] lists the ratios between the two main lines for each 
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Figure 4. Spectra of each maser listed in Table [4] arranged in order of decreasing R.A. Each maser is seen as two or more adjacent 
channels with a flux greater than three times the rms noise level. The data are smoothed with a Gaussian of width 2 channels. The x-axis 
velocity scale is calculated relative to the rest frequency of the 1665 MHz line. Continued on next page. 



maser from the observations discussed here, with limits given 
(based on the 3<r detection limit) where a detection was 
made in only one of the two main lines. 

In the MERLIN 1995 observations, four of the six 
masers are detected in both main lines, while 53.11+48.0 is 
detected only at 1665 MHz and 49.71+44.4 is present only 
at 1667 MHz. In 1997, the two masers detected were seen 
only at 1667 MHz. In 2002 six of ten are seen in both lines, 



53.11+48.0, 52.71+45.8 and 51.87+48.3 are seen only at 
1665 MHz, while 49.71+44.4 is again seen only at 1667 MHz. 

Of those which appear to be associated with Hn regions, 
most are brighter at 1667 MHz than at 1665 MHz. The only 
exception to this is 50.97+45.3 which is reversed only in 
the 1995 observations. The difference in magnitude between 
the two lines is, however, less than la. Two of the masers, 
53.11+48.0 and 52.71+45.8, appear to be associated with 
supernova remnants. These two masers are both detected 
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ID 


R.1995 


R2002 


Closest feature 


53.62+50.1 




1.02 


44.93+63.9 (Hn) 


53.11+48.0 


<0.47 


<0.41 


44.43+61.8 (SNR) 


52.71+45.8 




<0.15 


44.01+59.6 (SNR) 


51.87+48.3 




<0.45 




51.23+44.5 




1.22 


42.48+58.4 (Hn) 


50.97+45.3 


0.94 


1.42 


42.21+59.2 (Hn) 


50.51+45.6 




0.69 




50.36+44.2 


4.27 


4.72 


41.64+57.9 (Hn) 


50.02+45.8 


1.18 




41.29+59.7 (SNR) 


49.71+44.4 


>2.03 


>3.04 


40.96+57.9 (Hn) 


48.45+42.1 


1.67 


2.26 


39.68+55.6 (Hn) 



Table 5. Ratios of the measured intensities in the 1667 MHz line 
to the 1665 MHz line for each maser in M82. Limits are given using 
the measured noise level in each cube. The listed associations are 
the nearest co ntinuum features ( withi n a beam width of the maser 
position) from lMcDonald et al.l l|2002h , labelled according to their 
B1950 coordinates. 

only at 1665 MHz, if they are detected at all, with no emis- 
sion apparent at 1667 MHz. Further discussion of apparent 
associations with other known features are given in the next 
section. 

Investigating whether the masers vary in intensity be- 
tween the epochs is not as simple as comparing their flux 
measurements. The different channel widths used in each 
observation will determine, to some extent, the flux densi- 



ties measured. An alternative measure is the ratio between 
the two main line intensities in which the difference in chan- 
nel width will have the same effect on each of the two lines. 

It can be seen from Table [S] that, to within the uncer- 
tainties of the flux measurements, all of the ratios are consis- 
tent over the different epochs of observation. The limits on 
the ratio for 49.71+44.4 also vary by an amount larger than 
the uncertainties in the measured flux densities, with the 
1667 MHz line becoming comparatively bright in the 2002 
data. This source is interesting as it has additional struc- 
ture which is apparent in the 2002 observations (see Section 
153) . 

One source, 50.02+45.8, was detected in 1995, but is 
not visible in the 2002 observations. This source is one of 
the weakest seen in 1995 and may have faded below the 
noise level in subsequent observations. Several sources ap- 
pear in the 2002 data-set, but not in the 1995 data. One 
likely cause of this is due to the width of the channels in the 
1995 data-set which will average narrow emission so that 
it appears weaker, possibly reducing it below the detection 
thre shold. Variability is often observed in Galactic masers 
(e.g. ICaswell et al.1 1980), so it is possible that these masers 
are varying in intensity between observations. 

5.2 Associations 

Figure [5]shows the positions of both the m ain line maser de- 
tections described here and the SNRs from lMcDonald et al.l 
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Figure 5. The location of the main line masers (crosses) compared to the supernova remnants from McDonald ct al. (200^) (stars). The 
outline of the galaxy is provided by a 1— a contour of the continuum emission at 18 cm. 
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Figure 6. The location of the main line masers (crosses) compared to the Hn regions from iMcDonald et ail ll2002h (diamonds). The 
outline of the galaxy is provided by a 1— <j contour of the continuum emission at 18 cm. 
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Figure 7. The loc ation of the main lin e mase rs (crosses) compared to satellite line masers from [ Seaquist et al.l lll997T ) (broken crosses), 
water masers from IBaudrv fc Brouilletl l|l996h (squares) and the 2.2 /^m peak from [Lester et al .1 jl99Ch (triangle). The outline of the 
galaxy is provided by a 1— a contour of the continuum emission at 18 cm. 



(2002). Figure [6] shows the posit ions of the OH maser s 
along with the Hn regions from iMcDonald et all (2002). 
Figure [7] shows the positions of the main line OH masers 
de scribed here along wit h the satellite line masers detected 
bv lSeacmist et all l| 19971 ), the H 2 masers (both det e ctions 
and possible detections) lis ted inlBaudrv fc Brouilletl (|l996l ) 
and the 2.2 /im peak from iLester et all (|l990h . In all three 
of these figures, the outline of M82 is represented by a la 
contour from the continuum emission in the VLA 2002 ob- 
servation. From these figures is can be seen that most of 
the masers described here are closely associated with other 
known objects within M82. 

Four sources (51.23+44.5, 50.97+45.3, 50.37+44.1 and 
48.45+42.1) are all apparently coincident with a water maser 
feature to within the size of the beam in the VLA 2002 ob- 
servations. A further two (50.51+45.6 and 49.71+44.4) are 
located close to possible H2O maser detections. Of these 
associations, 50.51+45.6 is the most marginal as the separa- 
tion is of the order of the beam size. 51.23+44.5 is within a 
beam of two separate water m aser features, one of which is 
a marginal detection, listed in lBaudrv fc Brouilletl |l996h . 

Only a minority of the main line OH masers described 
here (four of eleven) do not appear coincident with water 
maser detections. Of these, all except 51.87+48.3 are coin- 
cident with other known fe atures. All of the de fi nite w ater 
maser detections listed by IBaudrv fc Brouilletl (|l996l ) are 
coincident with an OH main line maser detected in these 
observations. 

Six OH masers (53.62+50.1, 51.23+44.5, 50.97+45.3, 
50.36+44.2, 49.71+44.4 and 48.45+42.1) appear to be asso- 
ciated with Hn regions, while three (53.11+48.0, 52.71+45.8 



and 50.02+45.8) are within a beam of known supernova rem- 
nants. The large number of apparent associations is not sur- 
prising given the crowded fields and that maser emission is 
the result of an amplification of background continuum emis- 
sion. From the line ratios listed in the previous section it can 
be seen that the sources apparently associated with Hn re- 
gions are generally brighter at 1667 than at 1665 MHz, while 
those close to supernova remnants are brighter at 1665 MHz 
(apart from 50.02+45. 8). In studies o f Galactic masers such 
as that carried out bv lCaswelll l|l998h . most masers in star 
forming regions are found to be brighter at 1665 MHz than 
at 1667 MHz. This is discussed further in Section [6] 

The region around 49.71+44.4 contains two Hn regions 
as well as a water maser. The two Hn regions are located 
at R.A. = 49=707, Dec. = 45'.'03 and R.A. = 49:716, Dec. 
= 44'.' 13. Along with the multiple velocity components seen 
in this maser in the 2002 observations, this region of M82 is 
particularly interesting. Figure [8] shows this region in more 
detail. The location of this mas er also appears to be on the 
edge of one of the Hl/CO shells (|Wills et al.ll2002r i and there 
is a 1720 MHz maser detection within a beam width. 

The brightest X-ray source in M82 is located at 50! 2, 
46'.'7 (jKaaret et al.ll200ll ) with a la error circle of 0.7 arc- 
seconds. None of the maser detections lie within this circle, 
the closest detection being 50.02+45.8 which lies at a dis- 
tance of 2'.'85 (a linear distance of 44parsecs) so is clearly 
not associated. 

The black hole candidate described bv lPtak fc Grif fiths 
(1999) is located at R.A. = 09 h 55 m 52 s , Dec. = +69°40'48" 
with a 1.3 arcminute diameter error circle. The nearest 
maser detection to this source is 51.87+48.3, a source only 
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Figure 8. The region around 49.71+44.4 showing the OH ab- 
sorption from the VLA 2002 observation (dotted lines), the Hi 
shell (solid circle), Hn regions (diamonds), water maser (square), 
supernovae (stars), 1720 MHz maser (broken plus) and the OH 
maser position from the 2002 data-set (cross). The symbols are 
1" in size, apart from that for 49.71+44.4 which is the size of the 
restoring beam. 



present in the 2002 data-set at 1665 MHz. The separation 
between these positions is l'.'9 so considering the uncertain- 
ties it is possible that they could be ass ociated. 

The AGN can d idate suggested by IWills et all (|l997l ) 
and ISeaauist et alJ (|l997l ) is spatially coincident with the 
maser 52.71+45.8. This maser is apparently associated with 
an SNR, and one of only three detected at 1665 MHz while 
not detected at all at 1667 MHz. It is also the only one co- 
incident with a satellite maser from lSeaquist et al.l (|l997l ). a 
position which displayed absorption at 1720 MHz, but strong 
emission at 1612 MHz. 

For completeness, Figure shows the positions of the 
masers listed in Table [7] along with the po sitions of the Hn 
region s and supernova remnants list ed in McDonald et ail 
(2002), the satellite mas ers listed in ISeaauist et al. 

I (Il997t> . 

the Hi shells listed in IWills et alJ (|2002h . H?Q masers 
from [ Baudry fc BrouilletJ (|l996T ). and the 2.2 /im peak from 
iLester et alJ I 



5.3 Comparison with other maser detections 

The two brightest masers, 50.97+45.3 an d 50.36+44.2, are 
coinci dent with the OH masers reported bv lWeliachew et al.l 
l|l984h . corresponding to their masers m2 and ml respec- 
tively. In their data, one source (labelled ml) is seen at 
1665 MHz, while both their ml and m2 are detected at 
1667 MHz. In the 2002 VLA data, both masers are detected 
in both main lines at velocities of 80 and 98 ± 18 kms -1 re- 
sp ectively. This is consisten t with the velocities of the masers 
in IWeliachew et al' (1984) where their ml is at 65 to 80 
kms -1 , and m2 has a velocity of 85 to 100 kms - . 

Two of the OH masers reported here are appar- 
ently coincident w i th th e H2O features reported by 
iBaudrv fc BrouilletJ (|l996l ) to within the size of the l'.'4 



Table 6. List of the Gaussian components fitted to 50.97+45.3 
and 50.36+44.2 for channel 47 in the VLA 2002 data. Each com- 
ponent had the same major and minor axis, and position angle 
as the restoring beam (1.34x1.31 at 23°. 19). 



VLA beam in their observations. Their S2 is f/.'29 from 
the MERLIN position of 50.97+45.3 while the brightest 
water maser, S3, lies 0'.'63 from the MERLIN p ositio n of 
50.36+44.2. SI and S4 from lBaudrv fc Brouilletl (|l996l ) are 
both >l'.'6 from the nearest MERLIN detection. One other 
OH maser, 49.71+44.4, is located 0'.'57 from a possible wa- 
ter maser detection listed by the same authors. The veloci- 
ties of the H2O masers are consistent (within the uncertain- 
ties) with those measured for the OH detections presented 
here, apart from perha ps 50. 36+44. 2/S3 alt h ough even here 
the difference is small: IBaudrv fc BrouilletJ (|l996l ) measure 
103.9 kms - for S3 where as 50.36+44.2 is found to have a 
velocity of 80+18 kms -1 from the 2002 VLA data. 

Only one OH maser is coincident wit h a sa tellite line 
maser position reported by ISeaauist et al.l (|l997l ) : the sep- 
aration of 52.71+45.8 with feature 4 from that paper (a 
feature seen in absorption at 1612 MHz but strong emis- 
sion at 1720 MHz) is only 0'.'57. Both observations were car- 
ried out using the VLA with a #hpbw ~ l"4. This position 
is also coincident with the SNR 44.01+59.6, one of three 
definite dete ctions in this sample coincident with an SNR. 
Feature 4 in ISeacmist et al. (1997) had a measured flux of 
31.6+1.7 mjy at 1612 MHz in 1996 November and a peak of 
only 1.98 mjy beam -1 at 1665 MHz in 2002 November. At 
1720 MHz there is no emission apparent, only an absorption 
feature is present at this position. 



5.4 Spatial structure 

Most of the masers are unresolved in the observations de- 
tailed here. This is not surprising given that the typical size 
of Galactic maser spots is much less than a parsec, consid- 
erably less than the resolution of the observations detailed 
here. One maser feature does show extension in the VLA 
data which resolves into two components when observed 
with MERLIN. 

The extension to the SW of 50.97+45.3 in the VLA 
2002 data-set can be seen in Figure 1101 There is also a 
marginally detected second spatial component to this maser 
in the MERLIN 1995 data at the same position angle as 
the elongation seen in the 2002 data (see Figure lll[) . al- 
though it is only present at just over 3a. The weak detec- 
tion suggests that there is structure in the MERLIN data- 
set on scales larger than those to which the interferometer 
is sensitive which is resolved out. Using the AIPS task jm- 
fit to fit Gaussian components to a region including masers 
50.97+45.3 and 50.36+44.2 in the VLA 2002 data produces 
a better quality fit when three components are used, rather 
than two. The positions of the resulting fitted Gaussian com- 
ponents are given in Table [6] 




Figure 9. The location of the main line masers in M82 with other features shown for comparison. The solid line is a la contour of the 
continuum emission from the VLA data cube. The full crosses mark the po sitions of the main line masers reported here. The crosses 
with gaps mark the 1612/17 20 MHz masers reported in lSeaquist et all lll997h . The diamonds and stars mark the Hn regions and SNRs 
from iMcDonald et alj d2002h respectively. 
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Figure 10. Three of the central masers in a continuum- 
subtracted single channel from the VLA 2002 data-set showing 
the extension of maser 50.97+45.3. The position of the bright- 
est continuum feature (the SNR 41.9+57.5) is marked by a star. 
Contours are at (-1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10) X 1.751 mjy/beam. 



ID 


Velocity 


Nearest continuum 


(J2000) 


(kms" 1 ) 


feature (B1950) 


53.62+50.1 


347+9 


44.93+63.9 (Hn) 


53.11+48.0 


281+9 


44.43+61.8 (SNR) 


52.71+45.8 


233+9 


44.01+59.6 (SNR) 


51.87+48.3 


238+9 




51.23+44.5 


119+9 


42.48+58.4 (Hn) 


50.97+45.3 


96+9 


42.21+59.2 (Hn) 


50.51+45.6 


203+9 




50.36+44.2 


79+9 


41.64+57.9 (Hn) 


49.71+44.4 


148+9 


40.96+57.9 (Hn) 


48.45+42.1 


107+9 


39.68+55.6 (Hn) 



Table 7. The velocities of the ten definite maser detections in the 
VLA 2002 data. The features are named according to their J2000 
positions (relative to 09 h 55 m +69°40'). Velocities are measured 
with respect to the line rest velocity and where a maser is de- 
tected in both lines the velocities are consistent between the two. 
The nearest continuum features are from, and are label led ac- 
cording to, the convention used in iMcDonald et al.l (I2002T) (their 
corresponding J2000 IDs are given in Table [5] 



< 58.5 





Figure 11. One of the central masers in a single channel from 
the MERLIN 1995 data-set showing the extension of maser 
50.97+45.3. This image shows both the maser and continuum 
features. The position of the brightest continuum feature (the 
SNR 41.9+57.5) is toward the bottom of the image. Contours are 
at (-1, 1, 2, 4, 8, 16, 32 and 64) X 2.6mjy/beam. 



The maser 49.71+44.4 appears to consist of multiple 
components at different velocities in the 1667 MHz line. This 
velocity structure is seen in the MERLIN 1995 and VLA 
2002 data and may consist of multiple spatial components 
within a small region. This source is discussed in more detail 
in Section [531 



velocities of each feature measured in this data-set are given 
in Table [7] 

The velocities measured for most of the masers lie on a 
line consistent with the rotation of a disk; this can be seen in 
Figured! The sources 53.62+50.1, 53.11+48.0, 52.71+45.8, 
51.87+48.3 and 48.45+42.1 are at the same velocity as 
the OH absorption from the galaxy. On the other hand, 
51.23+44.5, 50.97+45.3, 50.36+44.2 and 49.71+44.4 lie on 
an arc which is blue shifted with respect to the absorption at 
the same position. The location of this arc is coincident with 
the po ssible shell seen in CO emission b y Matsushita et all 
(2000) and Hi absorption by IWills et all (|2002h (see Fig- 
ure 1 13p and could be associated with a superbubble a round 
the unusual source 41.95+57.5 l|Muxlow et al.ll2005l ). Fig- 
ure [3] shows that there is little OH absorption around this 
position. Another possibility is that these features, rather 
than being associated with a bubble or shell-like feature, 
are situated on the line o f the inner orbits of the possible 
bar structure described by I Wills et all (|2000h . 

One of the masers associated with this velocity feature, 
49.71+44.4 is unusual. In the VLA 2002 data-set there is an 
emission feature at 2.34 mjy in the 1667 MHz line with a ve- 
locity of 151 kms -1 . This feature has a weaker counterpart 
in the 1665 MHz line with a peak of <1 mjy at a velocity of 
168 kms -1 . There also appears to be a possible very weak 
second emission feature in the 1667 MHz line with a velocity 
of 222 kms -1 . This second feature is comparable with the 
noise in the map so was initially ignored. However, this sec- 
ond component in the 1667 MHz line is also present in the 
MERLIN 1995 data. 



5.5 Velocity distribution 

Figure [4] shows the spectra of the ten definite masers de- 
tected at greater than 5a in the 2002 observations. Most of 
the masers are clearly visible above the absorption but it 
can be seen that others, like 49.71+44.4, may be hidden in 
the absorption. It can also be seen that while some of the 
masers have the same line-of-sight velocity as the underly- 
ing absorption, several are blue or red shifted slightly. The 



6 COMPARISON WITH OTHER SURVEYS 

iBall fc Staehiil { 1968) compare OH maser detections at var- 
ious frequencies with other objects. They find that Hn re- 
gions are generally brig htest in either of the main lines (class 
I in the terminology of lTurne"rll 19691 ) while those associated 
with SNRs or other non-thermal sources associated with Hn 
regions are brightest at 1612 or 1720 MHz (Turner's class II) 
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Figure 12. A comparison between the maser velocities and the absorption. Crosses represent measurements of maser velocities from 
the 2002 data. The contours represent the OH absorption (and emission) from the VLA 2002 dataset and the horizontal dashed line 
represents the systemic velocity of M82. 



and suggest that the processes in non-thermal sources could 
inhibit the main line emission. 

In their conclusions, IPavlakis fc Kvlafisl |l996) state 
that where both 1665 and 1667 MHz lines are spatially co- 
incident, an FIR radiation field (assumed to be due to ther- 
mal emission from dust) must be present. They also state 
that the 1665 MHz line is typically stronger than that at 
1667 MHz. Most of the masers seen M82 appear to display 
the opposite behaviour, with the 1665 emission weaker, in 
most cases, than that at 1667 MHz. 

A survey bv lCaswelil (Il998l ) of OH ground state masers 
in the Southern Galactic plane confirmed previous conclu- 
sions that the 1665 and 1667 MHz lines often occur to- 
gether in space. The 1665 MHz line is usually stronger than 
the 1667 MHz line by a factor of around 3, although the 
1667 MHz line is stronger in some sources by up to a factor 
of 4. 

ICaswell et all (jl98Ch survey a selection of the Galactic 
plane for OH in a restricted longitude range. They detect 
a total of 40 main line maser sources, most of which are 
brighter at 1665 than 1667 MHz. Most are also associated 
with known Hn regions, while two of the eight remaining 
sources with no known Hn associations may be associated 



with SNRs. Three of the masers in M82 are similarly as- 
sociated with SNRs rather than Hn regions, although the 
possibility that there are undetected Hn regions along the 
same line of sight cannot be discounted. 

IWilson fc Barrettl (Il972l ) (and later ISlootmaker et all 
Il985h searched for OH masers associated with late- type stars 
in the Milky Way and found that most sources displaying 
emission in the main lines were brighter at 1667 MHz than at 
1665 MHz. They find that a positive 1667:1665 brightness ra- 
tio appears to be more characteristic of evolved stars, rather 
than star formation and Hn regions, with the differences in 
strength of emission due to the pro perties of the dust grains 
in the masing cloud (lElitzurll 19921 ). 

The most surprising result to come from this study is 
that, contrary to the results of other surveys discussed here, 
most of the masers in M82 are brighter at 1667 MHz and 
appear associated with Hn regions. This would imply that 
most of these masers are associated with regions of evolved 
stars, rather than Hn/OH regions as they appear to be. 
However, the apparent associations here are not necessarily 
physical and could just be due to line of sight coincidences 
since maser emission is the result of the amplification of 
background continuum emission. Another possibility is that, 
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Figure 13. The positions of the masers in the centre of M82 in relation to the Hi absorption shell seen bv lWills et al] {2002) at a velocity 
of 173 kms" 1 . Contours arc of the Hi absorption at (-10 to +10 in steps of 1) X 1.1 mjy/beam. 



as the resolution of these observations is 3pc at best, they 
are not physically associated with the observed Hn regions 
but with nearby regions of evolved stars. Higher angular 
ressolution measurements would be required in order to in- 
vestigate whether the associations are real and physical, or 
just due to line-of-sight effects. 



7 CONCLUSIONS 

Nine new OH main line masers have been detected in the 
central starburst region of the nearby galaxy M82. On a 
position-velocity plot of the galaxy, most of the masers lie 
along the same axis as the gas. Four of the masers are blue- 
shifted from the main distribution and could possibly be on 
the edge of an expanding shell or caused by gas orbiting 
within a bar. Most appear associated with Hn regions but 
are more characteristic of masers associated with late-type 
stars. 

Due to the spectral resolution of the 2002 VLA observa- 
tions the masers were not resolved in frequency, with most 
only visible in two channels. As the individual spots are 
an order of magnitude more luminous than typical Galactic 
masers, it is likely that they are made up of more than one 
individual masing region, so would have structure both spa- 
tially and in frequency. In order to investigate this further, 
observations have recently been made using the European 
VLBI Network in order to provide higher spatial resolution 
(much less than 1 pc at the distance of M82) , and the VLA 
at higher spectral resolution. It is hoped that further inves- 



tigation using these two observations will show evidence of 
both spatial and velocity structure. 

It is likely that there are more main line OH masers in 
M82 but, due to the depth of absorption and the low velocity 
resolution, faint or narrow masers could be buried to the 
extent that they are undetectable in the 2002 observations. 
The recent high spectral resolution VLA observations will 
have a gain in sensitivity over the 2002 observations which 
should enable the detection of weaker masers. 

Most of the maser spots detected are unresolved and 
likely consist of large numbers of small masing clouds. Any 
variability could be due either to variability of a few strong 
sources in a particular spot, the creation/destruction of 
many weaker maser sources within a spot, or a combina- 
tion of both. 

The association of most maser features with Hn regions 
is not surprising given the results of surveys carried out in 
out own Galaxy, although the comparative strength of the 
1667 MHz line is more characteristic of late-type stars than 
Hn regions. The association of two masers with known su- 
pernova remnants, while unusual, is not unique. A more sen- 
sitive survey for emission at 1612 and 1720 MHz would also 
help determine the nature of the maser sources. 
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ID 


1995 


1997 


2002 


Nearest feature 


Nearest feature 


(J2000) 


Sl665 


Sl667 


Sl665 


Sl667 


Sl665 


Sl667 


(B1950) 


(J2000) 


53.62+50.1 


<2.6 


<2.6 


<17 


<18 


2.44 


2.49 


44.93+63.9 (Hn) 


53.65+50.0 


53.11+48.0 


5.57 


<2.6 


<17 


<18 


1.86 


<0.77 


44.43+61.8 (SNR) 


53.14+47.8 


52.71+45.8 


<2.6 


<2.6 


<17 


<18 


5.02 


<0.77 


44.01+59.6 (SNR) 


52.73+45.7 


51.87+48.3 


<2.6 


<2.6 


<17 


<18 


1.71 


<0.77 






51.23+44.5 


<2.6 


<2.6 


<vr 


<18 


4.85 


5.93 


42.48+58.4 (Hn) 


51.22+44.5 


50.97+45.3 


8.35 


7.83 


<17 


23.0 


13.0 


18.5 


42.21+59.2 (Hn) 


50.95+45.2 


50.51+45.6 


<2.6 


<2.6 


<17 


<18 


1.85 


1.28 






50.36+44.2 


7.01 


29.9 


<17 


57.2 


10.7 


50.5 


41.64+57.9 (Hn) 


50.39+44.1 


50.02+45.8 


4.30 


5.09 


<17 


<18 


<0.77 


<0.77 


41.29+59.7 (SNR) 


50.06+45.9 


49.71+44.4 


<2.6 


5.29 


<17 


<18 


<0.77 


2.34 


40.96+57.9? (Hn) 


49.71+44.1 


48.45+42.1 


4.15 


6.94 


<17 


<18 


2.64 


5.97 


39.68+55.6 (Hn) 


48.44+41.9 



Table 8. All masers detected in the data-sets discussed here. Fluxes are jriven for each masers in each line in which is was detected in 



each epoch, along with the closest continuum feature from McDonal d et all ll2002f ) given in both B1950 and J2000 coordinates 
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